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(Dated: 13 August 2003) 

Abstract 

We present the results of a search for the production of new particles decaying into two jets in 
pp collisions at -y/i = 1.8 TeV, using the D0 1992-1995 data set corresponding to 109 pb~^. We 
exclude at the 95% confidence level the production of excited quarks {q*) with masses below 775 
GeV/c^, the most restictive limit to date. Wc also exclude standard- model-like W' {Z') bosons 
with masses between 300 and 800 GeV/c^ (400 and 640 GeV/c^). A W' boson with mass < 300 
GeV/c^ has been excluded by previous measurements, and our lower limit is therefore the most 
stringent to date. 

PACS numbers: PACS numbers: 13.85.Rm, 14.80.-j, 14.70.Pw 
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The direct production of hadronic jets is the dominant contribution to high transverse 
momentum (px) processes in antiproton-proton [pp) coUisions. There are many extensions 
of the standard model that predict the existence of new massive objects (e.g., excited quarks, 
W and Z' bosons that couple to quarks and/or gluons and may be observed as reso- 

nant structures in the two-jet mass spectrum. The previous observation of W and Z bosons 
decaying into two jets in the UA2 experiment jsf proved the feasibility of doing dijet mass 
spectroscopy at pp colliders. Subsequently, the UA2 0| and CDF experiments searched 
for new resonances in the dijet mass spectrum, and set limits on their production within the 
context of different theoretical models. This paper reports on a search for such resonances 
in the two-jet mass spectrum ^j, Q] using the data collected at a center-of-mass energy of 
1.8 TeV with the D0 detector in 1992-1995, corresponding to an integrated luminosity of 
109 pb-^ 

Jet detection in the D0 detector primarily utilizes the uranium/liquid-argon calorime- 
ters that cover the pseudorapidity region | ?7 1< 4, where rj = — ln[tan(6'/2)] and 9 is the polar 
angle with respect to the proton beam. Jets are reconstructed using an iterative jet cone 
algorithm with a cone radius of TZ=0.7 in 7^-0 space |6], where is the azimuthal angle. 
Background jets from isolated noisy calorimeter cells and accelerator losses are minimized 
via jet-quality criteria |.6]. The transverse energy of each jet is then corrected P| for offsets 
due to the underlying event, noise, multiple interactions and pileup, the fraction of parti- 
cle energy showering outside of the jet cone, and calorimeter energy response to incident 
hadrons. 

For each event that passes the quality criteria, the inclusive dijet mass can be cal- 
culated, assuming that the particles within the jets are massless, using the relationship 
= 2E^E^[cosh.{Ari) — cos(A0)], where and E^ are the transverse energies of the two 
highest-i?r jets. Each event is then weighted by the inverse of the efficiency of the quality 
criteria. The pseudorapidities of the two leading jets are selected to be | 771,2 |< 1-0 and 
At] =\ r]i — i]2 \ < 1.6 in order to maximize the range of dijet mass at which the trigger is 
efficient. 

A single trigger was used to collect the 1992-1993 data, with an E^ threshold of 115 
GeV, for an integrated luminosity of 14.1 pb~^. During 1994-1995, the data were collected 
using four triggers, with uncorrected Et thresholds of 30, 50, 85 and 115 GeV, for integrated 
luminosities of 0.36, 4.8, 56.5, and 94.9 pb~^, respectively. After the jet-energy corrections. 
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these trigger samples are used to measure the dijet mass spectrum above mass thresholds of 
180, 250, 320, and 470 GeV, respectively, where each of the triggers is > 97% efficient. The 
resulting dijet mass spectrum is shown in Fig. ^ The widths of the mass bins are chosen 
such that all events in any bin are recorded using a single trigger, there were > 10 events 
per bin, and the bin width is approximately equal to the mass resolution. 

The uncertainty in the dijet mass spectrum from the uncertainty in luminosity is 5.8%, 
and the uncertainty from the jet-quality and vertex criteria is 1%. The uncertainties due to 
the jet energy scale ^ are 7% (30%) for the lowest (highest) mass bin, and are correlated. 
The uncertainty in energy scale has three components: the uncorrelated, fully correlated, 
and partially-correlated uncertainties. A correlation matrix is calculated for the partially 
correlated uncertainties using the method described in Ref. jfl]. The uncertainties in the 
mass spectra due to the jet energy resolution are (0.5-3.0)% over the mass range under 
consideration. 

We consider three models for a possible signal in the diiet mass spectrum. The first 
model contains a mass-degenerate family of excited quarks l| that decay to a quark and 
a gluon (g* — > qg). We assume that the coupling parameters of the excited quarks equal 
unity if = f = fs = 1) and that the compositeness scale equals the mass of the excited 
quark (A* = Mq*). The second and third models 0\ contain additional W and Z bosons, 
respectively, with standard-model-like couplings, where all possible quark decays are allowed 
{W qq', Z' qq, with W' —>■ tb, and Z' tt allowed when kinematically possible). 
The leadingj-order W and Z' boson production cross sections are corrected by NLO "ii' 
facto." m of approximately L3, to account fo. Inghe.-otde. effect. All ,nodel= wete 
generated using pythia 6.2 [11|, with the CTEq6 parton distribution functions (PDFs) 



121. 

For each of the models, a Monte Carlo mass spectrum was generated at 25 GeV/c^ 
intervals from a mass of 150 GeV/c^ to 1 TeV/c^. Jets are reconstructed at the particle 
level using the same iterative jet cone algorithm that is applied to the data. The resulting 
energies are then smeared with the measured jet resolutions. Each of the mass spectra 
contains 50,000 events. Examples of the spectra generated for a resonant mass of 500 
GeV/c^ are shown in Fig. |21 

The data were analyzed using a Bayesian technique, with a fiat prior for the sig- 
nal (see Ref. [l3|). The predicted number of events per bin is given by /ij = 
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FIG. 1: The inclusive dijet mass spectrum. The events from each trigger have been corrected by the 
trigger's luminosity and event efficiency. The data were collected using triggers with uncorrected 
Et thresholds of 30 (open circles), 50 (solid squares), 85 (open triangles), and 115 GeV (solid 
stars). The error bars represent statistical uncertainties. 

{(^QCD^CgcDi + NxiCTxCxi) i^i^i where ctqcd^ is the predicted QCD two-jet cross section 
for mass bin i; Nx^ is the fraction of signal events in the bin {J2^Xi — 1); ctx is the 
cross section for the signal; Ci is the integrated luminosity; corresponds to the prod- 
uct of the efficiencies of the jet-quality criteria, the vertex selection efficiencies, and the 
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FIG. 2: The line shapes of a 500 GeV/c^ q* , W and Z' bosons , smoothed and normahzed to 
unit area. 

trigger efficiencies per bin; and CqcDi and are the jet energy and resolution correc- 
tions on the background and signal, respectively. Assuming Ni follows Poisson statis- 
tics, the probability that Ni events are observed in a given mass bin is then given by 
P{Ni I (Tqcd,,(^x, Nx,, C,ei,Cx„CQCDa I) = e'^^'fn^' /Nil, where / reflects all other "nui- 
sance" parameters. The probability of observing the set Ni that makes up the mass spectrum 
is then given by the product of these probabilities. To calculate the probability distribution 
for ax, Bayes' theorem is applied with the following assumptions about the prior probability 
distributions: ax has a uniform prior; aQCD^, ^i, CQCDi, Cxi and Ci all have Gaussian priors 
with widths given by their uncertainties; and Nx^ has a Poisson prior. 

Multiiet background was simulated using the next-to-leading order (NLO) program JE- 



TRAD with the cteq6m [Ij] PDFs, and renormalization scale (/i) of O.^E^^"^, where 
the E^^^ is the Et of the highest-i?T parton. Partons within 1.37^ of one another are clus- 
tered into a single jet if they are within TZ = 0.7 of their ii^-r- weighted ri,(f) centroid j^. The 
two highest-i?T jets are used to calculate the dijet mass, which is then smeared using the 
measured mass resolutions. The resulting distribution is shown in Fig. |S1 

A comparison between the background prediction and the data is given in Fig. 0] (only 
uncorrelated uncertainties are shown). The of the comparison is 25.0 for 25 degrees of 
freedom. This fit shows no obvious evidence for the existence of new particles. 



9 




10 

200 400 600 800 1000 1200 
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FIG. 3: The Jetrad (solid line) simulation of the inclusive dijet mass spectrum. The dashed- 
dotted lines show Pythia simulations of the excited quark line shapes for Mq* = 300, 500, 700, 
and 900 GeV/c^. 

The 95% confidence level (CL) limits on the production cross sections for the three 
resonances are extracted using the same Bayesian method described above. In Fig. El we 
compare our measured 95% CL limits (stars) with the expected cross section multiplied by 
the branching fraction {B) and acceptance for particles decaying to dijets (dashed curve). 
Branching fractions to all possible quark and gluon states are included in the acceptance. 
The acceptances for excited quarks {W and Z' bosons) range from 20% at 200 GeV/c^ 
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FIG. 4: The difference between data and the smeared Jetrad NLO QCD prediction normaUzed 
to the theoretical prediction ((Data — Theory) /Theory) using the cteqGm PDFs and a single 
renormalization scale /i = 0.5E^^^. The vertical error bars represent the sum of the uncorrelated 
uncertainties added in quadrature, while the horizontal error bands represent the widths of the 
mass bins. The highest mass bin extends to 1400 GeV/c^. 

to 60% (50%) for masses above 700 GeV/c^. We exclude excited quarks with Mg* < 775 
GeV/c^. This is is the most restrictive limit on excited quark production to date. A W boson 
is ruled out in the mass range 300 < Mw' < 800 GeV/c^. Previous measurements I3| 
have excluded a W boson with mass below 300 GeV/c^; our new measurement therefore 
sets a far more stringent lower limit on a W boson mass of 800 GeV/c^. A Z' boson with 
mass between 400 and 640 GeV/c^ is also excluded. 
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FIG. 5: The 95% CL on the production cross section multiplied by 5(X — > dijet) and acceptance, 
using the CTEQ6M PDFs for: (a) an excited quark q* (stars), compared with the predicted cross 
section (dashed line); Mq* < 775 GeV/c^ is excluded; (b) similarly, for a W' boson (stars), 300 < 
Mw' < 800 GeV/c^ is excluded; and (c) for a Z' boson (stars), 400 < M^/ < 640 GeV/c^ is 
excluded. 
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